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Abstract 

We obtain a new exact solution to the field equations in the EGB modified theory of gravity for 
a 5-dimensional spherically symmetric static distribution. By using a transformation, the study 
is reduced to the analysis of a single second order nonlinear differential equation. In general the 
condition of pressure isotropy produces a first order differential equation which is an Abel equation 
of the second kind. An exact solution is found. The solution is examined for physical admissability. 
In particular a set of constants is found which ensures that a pressure-free hypersurface exists 
which defines the boundary of the distribution. Additionally the isotropic pressure and the energy 
density are shown to be positive within the radius of the sphere. The adiabatic sound speed 
criterion is also satisfied within the fluid ensuring a subluminal sound speed. Furthermore, the 
weak, strong and dominant conditions hold throughout the distribution. On setting the Gauss- 
Bonnet coupling to zero, an exact solution for 5-dimensional perfect fluids in the standard Einstein 
theory is obtained. Plots of the dynamical quantities for the Gauss-Bonnet and the Einstein case 
reveal that the pressure is unaffected while the the energy density increases under the influence of 
the Gauss-Bonnet term. 
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I. INTRODUCTION 


Alternate or extended theories of gravity have aroused considerable interest recently in 
view of difficulties with the general theory of relativity to explain anomalous behaviour of 
gravitational phenomena such as the late time expansion of the universe. The mathematical 
reason for this interest is that the higher order derivative curvature terms make a nonzero 
contribution to the dynamics. In particular, Einstein-Gauss-Bonnet (EGB) theory has 
proved promising in this regard, and therefore is the most extensively studied. It appears 


Q. 


in a natural way in the effective action of heterotic string theory in the low energy limit 
Several new results have been reported especially dealing with the aspect of gravitational 
collapse. The causal structure of the singularities is different from general relativity for 
inhomogeneous distribution of dust and null dust ^]. 

Historically black hole models in EGB theory have been intensively studied. Boulware and 
Deser jf]] generalized the higher dimensional solutions in Einstein theory due to Tangher- 
lini [jj], and by Myers and Perry jf| to include the contribution of the EGB theory with 

8] have 


quadratic curvature terms. Wheeler j(j] , Torii and Maeda 121 and Myers and Simons 


0 


also considered black hole solutions in EGB theory. Inhomogeneous collapse of dust, that 
is pressure-free fluid with non-interacting particles, was studied by Maeda 9]. However 
explicit exact solutions were obtained by Jhingan and Ghosh [ic|. Dadhich et al 11] proved 
that the constant density Schwarzschild interior solution is universal in the sense that it is 
valid in both higher dimensional Einstein theory as well as in EGB gravity. The matching of 
these exterior metrics to an interior for brane world stars to produce analytical models was 


investigated by Casadio and Ovalle 


12 ] 


exterior was analysed by Clifton et al 


13] 


'lie matching of isolated masses to a Schwarzschild 


To date, there appears to exist no complete stellar model in 5-dimensional EGB theory for 
the perfect fluid configuration of static spherically symmetric matter. The constant density 


configuration of Dadhich et al ]11] and the static spherically symmetric star of Kang et al 
jr3 ] need to satisfy the junction conditions of EGB gravity so that matching is possible at 
the stellar surface. The junction conditions for EGB were derived in Davis 15] which are 
nontrivial and very different from general relativity. Note that the variable density model 
of Kang et al [l4] requires a further integration to produce an exact solution. As far as we 
are aware there is no known interior variable density spherically symmetric exact solution to 
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the EGB field equations. This is the object of our study in this paper and to highlight the 


role of junction conditions of 


15]. The problem of finding an exact interior metric in EGB 


theory for spherically symmetric distribution reduces to solving a system of three nonlinear 
partial differential equations in four unkowns: the dynamical quantities pressure and energy 
density and two gravitational potentials. As the system is under-determined, it is necessary 
to specify an additional contraint in order to close the system. This process is analogous to 
the standard Einstein gravity case. Traditionally, the approach has been to specify one of 
the four unknowns and by integration to resolve the remaining three. The reader is referred 


to the comprehensive listing of exact solutions found in this way by Stephani et al m and 
Delgaty and Lake [l7] for the Einstein case. When one of the field equations is replaced by 
the equation of hydrodynamic equilibrium, namely the vanishing divergence of the energy 
momentum tensor, then it is prudent to invoke an equation of state relating the energy 
density and pressure. This appears to have been the approach of 14[in their attempt to 
find an exact model. Interestingly the interior model presented in 14] generates the well 
known vacuum metric Q] of EGB in the limit of vanishing pressure and density. 

Recently Izaurieta and Rodriguez 18] argued that four dimensional gravity may effec¬ 
tively emerge from 5-dimensional EGB theory. The addition of diffeomorphism-invariant 
terms to the action principle leads to second order equations of motion and are therefore 
physically palatable. Consequently investigations in five dimensions are laboratories for ex¬ 
amining the impact of extra dimensions on physics, for example by demonstrating that a 
new exact solution satisfies elementary physical properties demanded of astrophysical ob¬ 
jects. In four dimensions the higher order curvature terms in EGB theory do not affect 
gravity. It is only with spacetime dimensions five or greater when the Gauss-Bonnet term 
contributes nontrivially to the dynamics. The simplest case in higher dimensions is five 
which has been extensively studied in several physical scenarios. We point out that the 
addition of the extra spatial dimension has been investigated by Kang et al in static 
stars, Brihaye and Reidel 19] in rotating boson stars, Ghosh et al Q in spherical collapsing 
bodies, and Chervon et al 20] in emergent universe models. The presence of additional 
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dimensions may have a dramatic effect on the behavior of matter. For example Maeda 
showed that massive timelike naked singularities may exist in five dimensions and massless 
ingoing null naked singularities are formed in dimensions greater than five in EGB theory. 
Massive timelike naked singularities do not exist in general relativity. Also note that the 
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dynamics of charged radiating gravitational collapse of shear-free matter has been studied 


recently in the socalled modified Gauss-Bonnet gravity 


21 ]- 


Our intention in this paper is to solve the nonlinear EGB equations for a static spherically 
symmetric matter distribution. In section II we briefly outline the basic equations in EGB 
gravity. The field equations in 5-dimensional EGB gravity are presented for a spherically 
symmetric metric, and they are then transformed to an equivalent form in Section III. An 
exact solution to the EGB equations is found in Section IV. In Section V a corresponding 
exact solution, in the Einstein case, to the 5-dimensional case is presented. The physical 
features of the model are investigated in Section VI. Some concluding remarks are made in 
Section VII. In the Appendix we present several exact solutions in the Einstein case in five 
dimensions by specifying a form for one of the metric potentials. 


II. EINSTEIN-GAUSS-BONNET GRAVITY 


The Gauss-Bonnet action in five dimensions is written as 


S 


V^g 


- (R — 2A + aL GB ) 


d 5 x + S ma tten 


( 1 ) 


where a is the Gauss-Bonnet coupling constant. The strength of the action Lqb lies in the 
fact that despite the Lagrangian being quadratic in the Ricci tensor, Ricci scalar and the 
Riemann tensor, the equations of motion turn out to be second order quasilinear which are 
compatible with a theory of gravity. The Gauss-Bonnet term is of no consequence for n < 4 
but is generally nonzero for n > 4. 

The EGB field equations may be written as 


Gab + O.H a b — T a b, (2) 

with metric signature (—h + + +) where G a b is the Einstein tensor. The Lanczos tensor is 
given by 

H a b = 2 ( RRab — 2 RacRb ~ 2 R cd R ac bd + R C a e Rbcde) ~ ^dabLcB, (3) 

where the Lovelock term has the form 


r _ 7 d2 i T) -pabcd, a e> r>cd 

J-'GB — A + Kabcd-ti ~ . 


(4) 
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III. FIELD EQUATIONS 


The generic 5-dimensional line element for static spherically symmetric spacetimes is 
taken as 

ds 2 = — e 2u dt 2 + e 2X dr 2 + r 2 [d6 2 + sin 2 9d<f> + sin 2 6 sin 2 (f)di(i) , (5) 

where v(r) and A(r) are the gravitational potentials. We utilise a comoving fluid velocity of 
the form u a = and the matter field is that of a perfect fluid with energy momentum 

tensor T a b = ( p + p)u a Ub + pg a b- Accordingly the EGB field equations (J2]) reduce to 

p = ^ 3 [4a\' + re 2X — re 4A — r 2 e 2X X' — 4ae 2A A') , (6) 

p = ^ 3 (—?’e 4A + ( r 2 is ' + r + 4 au') e 2X — 3ai/') , (7) 


P = 


g4A^2 


e 2\ r 2 


e 2\, r 2 


-e 4A — 4av" + 12av'M — 4a (z/') 2 j 
^1 — r 2 v' A' + 2rz/ — 2rA' + r 2 (z/') 2 j 
40?^'A' + 4a (v') 2 + 4az/"^ . 


r 2 v" 


( 8 ) 


Note that the system ([6])-((8]) comprises three field equations in four unknowns which is 
similar to the standard Einstein case for spherically symmetric perfect fluids. Observe that 
the vacuum metric describing the gravitational field exterior to the 5-dimensional static 
perfect fluid may be described by the Boulware-Deser [s.] spacetime as 

dr 2 


ds 2 = —F(r)dt 2 + 


Hr] 


+ r 2 (d6 2 + sin 2 9d(p + sin 2 9 sin 2 (j)d '?/>) , 


(9) 


where 


F(r) = 1 + — 
4a 


1 + 


8 Ma 


In the above M is associated with the gravitational mass of the hypersphere. The exterior 
solution is not unique and neither is there a Birkhoff theorem analogous to the 4-dimensional 
gravity case. At least this metric involves branch cuts. Bogdanos et al 22] have analysed 
the 6-dimensional case in EGB and demonstrated the validity of Birkhoff’s theorem for this 
order. 

We invoke the transformation e 2u = y 2 (x), e~ 2X = Z(x) and x = Cr 2 (C being an 
arbitrary constant) which was utilised successfully by Durgapal and Banerji [23], Finch and 
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Skea 24] and Hansraj and Maharaj 25] to generate new exact solutions for neutral and 


charged isotropic spheres. For applications to charged anisotropic rela 
this transformation see the recent works of Mafa Takisa and Maharaj 


dvistic matter with 


al 


26] and Maharaj et 


27] in four dimensions. The held equations 


3Z + 


may be recast as 
3(Z-l)(l-4aCZ) 


x 


P_ 

C' 


( 10 ) 


3(Z-1) | 6 Zy 24aC(Z-l)Zy _ p_ , n , 
x y xy C 1 


2 xZ (4aC[Z — 1} — x)y — (x 2 Z + 4aC 


xZ-2Z + 2Z 2 - 3xZZ 


y 


1 + xZ — Z ) y = 0, 


( 12 ) 


where the last equation is called the equation of pressure isotropy. Equation (fl2|i has been 
arranged as a second order differential equation in y, which for some analyses in the 4— 
dimensional Einstein models, proves to be a useful form. Functional forms for Z(x) may 
be selected a priori so as to allow for the complete integration of the held equations. For 
example, the form Z = 1 + x produces a higher dimensional Schwarzschild solution with 
constant density. This corroborates the result of Dadhich et al [ll] that the constant density 
Schwarzschild solution is universal - that is it is independent of dimension. We have also 
found a number of other cases for Z for which (IT2]) is integrable and these will be dealt with 
in the future. 

For the present work it should be noted that (1T2T) may also be regarded as a hrst order 
ordinary differential equation in Z, and may be expressed in the form 


( x 2 y + xy + 4aCxij — 12 aCxyZ) Z + 8 aC ( y — xy) Z 2 
+ (2 x 2 y + 8 aCxy — 8 aCy — y) Z + y = 0. 


(13) 


This is an Abel equation of the second kind for which few solutions are known. However note 
that by choosing forms for y should in theory result in the expressions for Z by integration. 
Therefore we seek choices for the metric potential y which will allow for a complete resolution 
of the geometrical and dynamical variables. 
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IV. NEW EXACT INTERIOR SOLUTION IN THE EGB CASE 


Locating exact solutions for (fl3l) is difficult to achieve in view of its nonlinearity. One 
strategy is to investigate the consequence of one or more of the coefficients to vanish. In 
requiring that the coefficient of Z 2 vanishes we obtain the restriction 

y~xy = 0, (14) 


which may be solved to give 

y = \ c ^ 2 + c 2, ( 15 ) 

where C\ and are constants of integration. Note that the restriction (fTlj) simplifies (I13p 
but does not remove its nonlinearity. Inserting (TT5|) into (TT31) the condition of pressure 
isotropy is transformed to 


[3Cix 3 + 8f5C\x 2 + 2 C 2 x - 24 f3C lX 2 Z] Z + [3Cix 2 - 2 C 2 ] Z + C x x 2 + 2 C 2 = 0, (16) 

where for convenience we set j3 = aC. Renaming = e equation (fIBl) assumes the simpler 
form 

[3ex 3 + 8 (3ex 2 + 2x — 24 /3ex 2 Z] Z + [3ex 2 — 2\ Z + ex 2 + 2 = 0, (17) 


which will aid our graphical investigations. The parameters e and j3 will have to be assigned 
values to obtain the qualitative features of the eventual model. On solving (fl7l) . we obtain 
the solutions 


3ex 2 + 8f3ex + 2 ± T 
24 /3ex 


(18) 


where 


T = V3e 2 (3x 4 + 32/3x 3 ) + 4e (16^ 2 e + 3 + 144/? 2 eC' 3 ) x 2 - 4(16 pex - 1), 


and C 3 is an integration constant. 

With the help of (ll8ll and (ITU the energy density for the EGB case is given by 
p 27e 2 x 4 — 48/3e 2 x 3 — 32f3 2 ex + 4 — 4(1 — 4 (3ex) T 
C 48(3e 2 x 4 

_x(3ex 2 + 16/3ex - 2) T' - (T - xV) T 

48 f3e 2 x 4 

while the pressure p has the form 

p 27e 2 x 4 + 96 e 2 /3x 3 + 8e (3 + 32 e/3 2 ) x 2 — 64 e/3x + 4 
C 24 e/3x 2 (ex 2 + 2) 

(3ex 2 + 16e/3x - 2 + 2T) T 
24e(3x(ex 2 + 2) 


(19) 
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( 20 ) 








via (ITT)) . The adiabatic sound speed parameter is found to be 


dp = 2 ex 2 U(x) 

dp (ex 2 + 2) 2 V(x) ’ 

where U(x) and V(x) are, respectively, given by 

U(x) = — 96e 3 /?x 5 + 4e 2 (15 — 128e/3 2 ) x 4 + 384e 2 /3x 3 — 16ex 2 
+ 128e/3x -16 + 8 (ex 2 + l) T 2 - 4x (ex 2 + 2) TT' 

+2 (3e 2 x 4 + 24 e 2 /3x 3 — 4ex 2 + 16 e/3x — 4) T 
—x (3e 2 x 4 + 16e 2 /3x 3 + 4ex 2 + 32/3ex — 4) T', 

and 

V(x) = 48e 2 x 3 + 96,5ex - 16 + 16 (1 - 3/3ex) T - 4T 2 
+x (3ex 3 + 48e/3 - 10) T' + x (5T + xT') T' 

—x 2 (3ex 2 + 16 e/3x — 2 + T) T". 

To study the energy conditions we need to obtain forms for p — p, p + p and p + 3 p. Explicit 
forms for these expressions can be immediately generated from (fl9l) and (j20jl . We will 
compare these expressions with the corresponding forms of the 5-dimensional Einstein case 
later. 

Other restrictions in (fT3l) may lead to new models in addition to that considered in this 
section. For example on setting the coefficient of Z in (fl3|l to zero, we obtain 

x(x + l)y-y-^y = 0, (22) 

where we have set 4 aC = 1 for simplicity. This produces a differential equation which is of 
the hypergeometric type, and consequently is not readily expressible in terms of elementary 
functions. Therefore the prospects of establishing Z explicitly are remote for the resultant 
form for y. This case may be treated with the other methods and this will be considered 
later. Additionally, it should be noted that the vanishing of the coefficient of Z is not 
mathematically feasible as it involves both y and Z. 



V. NEW EXACT INTERIOR SOLUTION IN THE EINSTEIN CASE 


The 5-dimensional Einstein version of the above is obtained by setting a = 0 in (fi~3l) to 
give the differential equation 


2 x 2 Zy + x 2 Zy + (1 — Z + Zx)y = 0. 


(23) 


This differs from its 4-dimensional counterpart 


4 x 2 Zy + 2 x 2 Zy + (1 — Z + Zx)y = 0, 


(24) 


in the coefficients of its first two terms 
reduces to 


23j. Invoking the metric ansatz (TT5]h equation 
x(3ex 2 + 1 )Z + (3ex 2 — 1 )Z + (ex 2 + 1) = 0, (25) 


where we have set e = The general solution to (1251) is given by 


2 c 2 ' 


Z(x) = 


1 + C 3 x — ex 2 


(26) 


1 + 3ex 2 

It is important to observe that this does not follow as a special case of solution (fl8|i since /3 
and consequently a appears in the denominator. Now the dynamical quantities, pressure, 
energy density, sound speed index, and energy conditions have the forms 
3 C [C 3 + 5C 3 ex 2 — 8e 2 x 3 ] 


P = 


P = 


(ex 2 + 1)(1 + 3ex 2 ) 
6 C [C 3 — 6ex — 6e 2 x 3 ] 


(1 + 3ex 2 ) 2 

dp x(l + 3ex 2 ) [12e 3 x — 13C 3 e 2 x A — 16e 2 x 3 — 6C 3 ex 2 — I2ex + C 3 ] 
dp 6(ex 2 + l) 2 [3e 2 x 4 + 6ex 2 — 2 C 3 x — 1] 


(27) 

(28) 
(29) 


P~P = 
P + P = 

p T 3p = 


3 C \C 3 - 12ex - 6C 3 ex 2 - 16e 2 x 3 - 15C 3 e 2 x 4 + 12e 3 x 5 ] 

(1 + ex 2 )(l + 3ex 2 ) 2 

3 C [3 C 3 - 12ex + 6C 3 ex 2 - 32 i 2 x 3 + 15C 3 e 2 x 4 - 36e 3 x 5 ] 
(1 + ex 2 )(l + 3ex 2 ) 2 

3C [5C 3 - 12ex + 26C' 3 ex 2 - 48e 2 x 3 + 45C' 3 e 2 x 4 - 84e 3 x 5 ] 


(30) 

(31) 

(32) 


(1 + ex 2 )(l + 3ex 2 ) 2 

respectively. In the plots to follow, we exhibit the above quantities in comparison with their 
EGB counterparts to investigate the role of the Gauss-Bonnet term in the solution. 
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VI. PHYSICAL FEATURES 


We extrapolate the familiar conditions imposed on stellar configurations in the usual 
Einstein theory to the 5-dimensional EGB case. It would be interesting to see if the models 
generated in this modified theory of gravity satisfy these standard conditions. We require 
that the energy density p and pressure p , and the metric potentials e 2 " and e 2A , should 
be regular in the interior. The radial pressure should vanish at the boundary r = R. The 
gradients p' and p' should be negative for barotropic matter. The speed of sound should 
remain subluminal throughout the interior of the star. At the boundary r = R the metric 
functions should match smoothly to the exterior Boulware-Deser jf| solution 


e MR) 


K 2 

4a 



= e -2\(R) 


For realistic matter we require compliance with the energy conditions: weak energy condition 
(p—p > 0), strong energy condition (p+p > 0), and dominant energy condition (p+Sp > 0). 

Clearly a complete analytic treatment of our solution is ostensibly not possible given its 
complexity. From an examination of the pressure (1201) it is evident that solving for x in 
terms of p is not attainable as this reduces to solving an eighth degree polynomial equation, 
and there is no general way of doing this. This in turn means that it is not possible to 
write the density as a function of pressure to obtain a barotropic equation of state; the form 
p = p(p) is not possible. 



FIG. 1: Plot of energy density versus radial coordinate x. 
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FIG. 2: Plot of pressure versus radial coordinate x. 



FIG. 3: Plot of sound-speed parameter versus radial coordinate x. 



FIG. 4: Plot of energy conditions versus radial coordinate x. 

We proceed to select a set of values for the various constants in the problem in order to 
determine a model that harmonises best with the physical conditions. The graphs in Fig. 
1-4 displayed were produced by assuming the values C\ = —0.02, C 2 = —50, C 3 = —0.05, 
a = —150, C = 0.0002 and (5 = aC. (Observe that the use of a negative coupling constant 
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a is not novel - see for example Guo and Schwarz 30].) For these choices it is pleasing 
to note from Fig. 2 that a pressure-free hypersurface does indeed exist and is defined 
approximately by x = 40.8647 geometric units. Within this boundary we note that the 
pressure is positive definite, and decreases monotonically outwards from the central regions 
to the boundary r = R. Additionally, Fig. 1 demonstrates that the density is positive and 
decreasing everywhere within the spherical distribution. Importantly we observe that the 
causality criterion 0 < ^ < 1 is satisfied everywhere in the interior of the star as evidenced 
by Fig. 3. This implies that the sound speed is never super luminal within the boundary. The 
energy conditions are depicted in Fig. 4. From these we infer that all the conditions: weak 
(solid line), strong (dashed line) and dominant (dotted line) are satisfied within the radius 
of the distribution. Using the values of the aforesaid constants as well as the boundary value 
x = 40, 8647 allows us to compute the mass as M = 1.51 x 10' geometric units. Working in 
these same units we find that the mass to radius ratio is -^ = 369700 which clearly violates 


the Buchdahl 


31] limit ^ < | valid for stars in the Einstein general theory of relativity. 
This suggests that the Buchdahl upper bound may not hold in this model in EGB gravity 
theory. 

The dashed curve in each of Fig. 1 and Fig. 2 reflects the situation when the EGB 
coupling constant is set to 0. In other words, it shows the 5-dimensional Einstein analogue. 
It can be seen from Fig. 2 that there is no discernible difference between the pressure profiles 
in the Einstein and the EGB scenarios. The density plot in Fig. 1 however, demonstrates 
that the radius for a positive energy density is improved by the presence of the EGB coupling 
constant. This implies that the gravitational held in the Gauss-Bonnet theory can sustain 
a greater amount of matter per unit radius as opposed to its Einstein counterpart. Note 
that the central singularity in this model is an artefact of the non-removable curvature 
singularity. A way of avoiding the singularity is to call upon a two-fluid scenario such as in 
a core-envelope model. This idea is popularly invoked in the standard 4-dimensional theory, 


see for example the deconhned quar 


core model surrounded by an envelope of barotopic 


matter of Sharma and Mukherjee 32], 


Matching of the interior metric with an exterior metric such as the Boulware-Deser solu- 
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tion is achieved via setting 


\c 1 c 2 r' + c^ - i + + 


8 Ma\ 


R 2 ( 8 Ma 

1 H-1 - i H- -r 

4a V R 4 


3C\C 2 R 4 + 8/3C'iC',R 4 + 2C 2 ± C 2 T 
2Af3C\CR 4 : 


(33) 

(34) 


0 = 27 C 2 C 4 R 8 + 96 C 2 PC 3 R 6 + 8Ci (3C 2 + 32Ci/5 2 ) C 2 i? 4 
-64C' 1 C 2 /3C' j R 2 + 4C*2 

- (3C'iC' 2 C' 2 i? 4 + \<oC x C 2 $CR 2 - 2 C 2 + 2C' 2 T) T. (35) 

The last equation arises as a result of the vanishing of the pressure at r = R. There are 
sufficient free parameters in the model to ensure that conditions (133|) (135|) are satisfied. We 
find that 


_ —b ± 2aE 1 
Cl “ 4a ± 2 


4 a 3 EB ± [4a6e - 6 3 - 8a 2 /] 


4 a 3 E 


(36) 


2 A - C^R 4 


(A + C X C 2 R 4 + AfiC^R 2 [1 - 37L 2 ]) 2 GA/AC^CR 2 - 4C 2 
27(3 2 CfC 2 R 4 (2A - C^R 4 ) 2 + A32f3 2 CfC 2 R 4 

9CjC 4 R 8 + 96 PC 2 C 3 R 6 + AC^R 4 (16/3 2 Ci + 3C' 2 ) 
A32/3 2 CjC 2 R 4 ’ 


(37) 


(38) 


where 


A 

a 

b 

e 


f 

h 


1 + 


IR 

Aa 


1-11 + 


8 Ma 


N 

HC 8 ^ 16 + 112fiC 7 R 14 + 2407l 2 /3C' 7 i? 14 - 192 ( 5 2 C' 6 i? 12 + 768 f3 2 C e R 12 , 


-196 AC 6 R 12 - 576 APC 5 R W - 1152A 3 /3C 5 R W - 10247l/3 2 C' 4 J R 8 - 3072 A 3 p 2 C 4 R 8 , 
—76A 2 C 4 R 8 - 16C 4 R 8 + 1088^ 2 C ,3 j R 6 + 17247l 4 C 3 i? 6 - 128 fiC 3 R 6 
+1792 A/3 2 C 2 R 4 - 256 f3 2 C 2 R 4 , 

68 AC 2 R 4 + 6A0A 3 CR: 2 - 128ACR 2 - A^fiCR 2 , 

64 A 4 , 
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with 


E 

B 


2^b 2 + ^3 x 2^b 2 — 2^ ae + Aag^j g 3 
\ 3 x 2 ia 2 g^ 

3 x 2§a (e 2 — 3 bf + 12a/i) + ^3 x 2^b 2 g^ — 2waeg3 — 2§a [e 2 — 36/ + 12a/i] j — 

3 x 2^a?g^ 


2 

2a^3 


and 


g = 2e 2 - 9bef + 27af 2 + 27b 2 h - 72aeh 

+ \/-A [e 2 - 3bf + 12a// 3 + [2e 2 - 9 bef + 27af 2 + 27b 2 h - 72aehf. 

From (BSD - dSHD we observe that the free parameters Ci, C 2 , C 3 in the model are defined in 
terms of a, M, R ( and C). It is interesting to note that the parameters that arise in the 
integration are defined in terms of physically relevant quantities: a is the Gauss-Bonnet 
coupling constant, R is the radius and M is the mass of the star. 


VII. CONCLUSION 


We have produced an exact solution for a static spherically symmetric distribution of 
perfect fluid in the modified EGB gravity theory. The model has been studied for physical 
admissability and has been found to satisfy several elementary tests for physical reality. 
The pressure function vanishes for a particular radius and this hypersurface identifies the 
boundary of the fluid. Within this boundary the pressure and energy density profiles are 
positive for the choice of constants made. Importantly the fluid is found to be causal as the 
sound speed is subluminal. The weak, strong and dominant energy conditions are found to 
be satisfied everywhere in the interior. Finally matching with the Boulware-Deser exterior 
metric is permitted. We point out that th e g eneral junction conditions for matching in 
EGB gravity has been considered by Davis 15]. The higher order curvature terms lead to 
a modified set of junction conditions to be satisfied at the stellar surface. These surface 
equations are very different from general relativity, and their complexity makes it unlikely 
to easily demonstrate an exact solution. For a complete stellar model of a star in EGB 
gravity those boundary equations should be satisfied. In this treatment we have shown that 
it is possible to find bound EGB interior solutions with variable densities and pressures with 
desirable physical features. 
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Appendix: Further New Exact Solutions in the Einstein Case 

We note that equation (fT2lh or m, admits large classes of solutions for prescribed forms 
of Z and y. A particular 5-dimensional exact model in the Einstein case (a = 0) was 
presented in Section IV. Other exact solutions to (TT3|) are possible for different forms of Z. 
In this appendix we list a few new solutions in various categories. However, for the sake of 
brevity we do not study these solutions in detail. 


1. The form Z = 1 + x r ' 


With this choice for Z, equation (1T21) assumes the form 


2x n (x n + 1 )y + nxy + (n — 1 )y = 0, 


(A.l) 


and has solutions which are the Legendre polynomials given by 


y(x) = Civ^LegendreP ( - 


-n 


+ Vn 2 - 12n + 12 1 


2 n 


,-,Vx n + l 

n 


(A.2) 


^ r - T , „ / -n+ y/ri 2 - 12 n + 12 1 , - 

+C 2 \/xLegendreQ ---, Vx n + 1 , (A.3) 

\ An n J 

where the standard Legendre functions of the first and second kind are defined, as usual, in 
terms of hypergeometric and gamma function. For certain values of n the Legendre functions 
reduce to elementary functions. Some of these cases are presented in Table 1. 

It should be noted that the n — 1 case corresponds to the model treated by Dadhich et 
al [ll], and which was shown to be equivalent to the interior Schwarzschild solution. The 
other cases appear to be novel. 
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TABLE I: Forms for the potential y for specific values of the parameter n. 


n 


Potential y 


1 


2C i \/l + x + C *2 


2 

- 1 
-2 

1 

2 



2. The form Z = (1 + x) 2 

The general form Z — (1 + x) n does not appear to yield elementary solutions except in 
the case n = 1 (coinciding with the above) and n = 2. In the latter case the exact solution 
is given by 


y(x ) = C\ cos 


ln(l + x ) 

WT 


+ C 2 sin 


ln(l + x) 


(A.4) 


where C\ and C 2 are integration constants. Consequently the line element for this solution 
has the form 

2 1 


ds 2 = 


C 1 cos 


ln(l + x) 

WT 


+ C 2 sin 


ln(l + x) 

WT 


dt 2 + 


(1 + x) : 


-dr 2 + r 2 dn 2 , 


where dQ 2 is the line element for the customary 3-sphere. 


3. The form Z = 


l+x 


: 5—dimensional Finch—Skea model 


Even though the case Z = is a simple form, it is worthy of special attention in its 
own right by virtue of the fact that this ansatz has been used by Finch and Skea 24] in four 
dimensions to generate physically reasonable stellar models that conform to the astrophysical 


theory of Walecka 


, and extended by Hansraj and Maharaj [25|] to include 


Finch-Skea model was a correction of the earlier work of Duorah and Ray 
5-dimensional case the master field equation (fl2|) is 


charge. The 


29 . For the 


2 (l + x)y-y + y = 0 . 


(A.5) 
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Note that its counterpart in 4-dimensions is given by 


4(1 + x)y - 2 y + y = 0, 


(A.6) 


where a difference in the initial coefficients is evident. The solution to (1A.5|) is given by 


y(x) = (C 2 + c\w) cos w + (c 2 w — ci) sinry, (A.7) 

where w = \/2(l + x ) and ci and C 2 are arbitrary constants to be established by matching 
with the exterior Schwarzschild-Tangherlini solution. On the other hand the Finch-Skea 
equation (IA.6[) is solved by 


y(x) = (c 2 + C\v) cosv + ( c 2 v — c\) sinu, (A.8) 

where v = \/l + x. Therefore it is clear that the physical properties of these solutions should 
be practically the same as there is only a difference of the factor of \/2. Clearly the change 
of dimension to five does not materially influence the physics of the perfect fluid. 
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